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Abstract 


Since the eariy 1990s the Aerofiightdynamics Directorate at the Ames Research Center 
has been conducting tests on isolated hingeless rotors in hover and forward flight. The 
primary objective is to generate a database on aeroelastic stability in trimmed flight for 
torsionally soft rotors at realistic tip speeds. The rotor test model has four soft inplane blades 
of NACA 0012 airfoil section with low torsional stiffness. The collective pitch and shaft tilt 
are set prior to each test run, and then the rotor is trimmed in the following sense: the 
longitudinal and lateral cyclic pitch controls are adjusted through a swashplate to minimize 
the 1/rev flapping moment at the 12% radial station. In hover, the database comprises 
lag regressive- mode damping with collective ^itch variations: —3° < 0 O < 8°. In forward 
flight the database comprises cyclic pitch controls, root flap moment and lag regressive- 
mode damping with advance ratio, shaft angle and collective pitch variations: 0 < fi < 0.36, 
0° < a s < 6°, and 0° < &o < 5.9°. This report presents the predictions and their correlation 
with the database. A modal analysis is used, in which nonrotating modes in flap bending, lag 
bending and torsion are computed from the measured blade mass and stiffness distributions. 
The airfoil aerodynamics is represented by the ONERA dynamic stall models of lift, drag and 
pitching moment, and the wake dynamics is represented by a state-space wake model. The 
trim analysis of finding the cyclic controls and the corresponding periodic responses is based 
on periodic shooting with damped Newton iteration; the Floquet transition matrix (FTM) 
comes out as a byproduct. The stability analysis of finding the frequencies and damping 
levels is based on the eigenvalue-eigenvector analysis ot the FTM. All the structural and 
aerodynamic states are included from modeling to trim analysis to stability analysis. A 
major finding is that dynamic wake dramatically improves the correlation for the lateral 
cyclic pitch control. Overall, the correlation is fairly good. 


N omenclat ure 

Unless otherwise stated, the symbols below are dimensionless: 
a linear lift curve slope 

a,-, o. m damping factors in dynamic stall drag and pitching moment models 



b 

airfoil semi-chord, (1/i?) 

c, c 

airfoil chord, (m, 1/i?) 

Cd, c do 

airfoil drag coefficient and constant profile-drag coefficient 

C d ■ 

quasisteady drag coefficient 

c z ,c m 

airfoil lift and pitching moment coefficients 

Cm „ 

airfoil pitching moment coefficient at zero angle of attack 

Cm, 

quasisteady pitching moment coefficient 

Ct 

thrust coefficient 

Ct/o's 

thrust level 

c« 

extrapolated linear-lift coefficient 

C z . 

quasisteady lift coefficient 

d 

dynamic-stall-lift damping parameter 

e 

phase shift parameter in dynamic stall lift model 

ey, e/ 

radial station for flapping moment estimation, (1 //?, m) 

e*, eh 

hinge offset, (1 / R, m) 

Ed, E m 

phase shift parameters in dynamic stall drag and 
pitching moment models 

FM}j , . . . , FA// 2 

Galerkin-type integrals, Eq. (20) 

I,J,K 

unit vectors associated with inertial frame 

i,j,k 

unit vectors associated with undeformed blade coordinate system 

' v mi ? rt 77i2 

mass radii of gyration of blade cross section about its principle axes 

/, i 

blade span, 1 — (1 //?, m) 

Ltd 

sectional aerodynamic forces 

c c 

total horizontal and vertical forces, (N) 

771 

mass per unit length, (kg/m) 

m r 

reference mass per unit length, (kg/m) 

M* 

aerodynamic pitching moment 

M v 

total pitching moment, (Nm) 

M 

total flap moment at e/, (1/ pbVt 2 R A ) 

Ad nc 

n-th harmonic cosine component of Ad 
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Mo 

r c 

rd,r m 


r i 

R 

t 

V 

V 

Vity) 

ID 


W 

Wity) 
X, Y, Z 

X, X 
a 


a* 


(%ss 

flpc 

7 

6 

0 

0o 

0a , 0s 

X 

A 2 

0 


n-th harmonic sine component of M. 
steady component of M. 
root cutout, (1 /R) 

frequency parameters in dynamic stall drag and 

pitching moment models 

radial station of the i-th blade, (1 / R) 

rotor radius, (m) 

time, ( 1 /n) 

inplane (lag) bending deflection 

dimensional inplane (lag) bending deflection 

generalized coordinate of i-th mode for lag degree of freedom 

dynamic stall lift frequency parameter; also out of plane 

(flap) bending deflection, (1 / R) 

dimensional out of plane (flap) bending deflection 

generalized coordinate of i-th mode for flap degree of freedom 

inertial coordinate system 

undeformed coordinate system 

radial distance measured from the rotor center, (1 / R, m) 

blade airfoil angle of attack 

shaft angle, (deg, rad) 

quasisteady stall angle, (deg) 

blade precone, (deg) 

Lock number (blade inertia parameter) 

pitch-rate coefficient 

blade pitch angle 

collective pitch angle, (deg, rad) 

lateral and longitudinal cyclic pitch angles, (deg, rad) 

time-delay parameter 

flap bending stiffness, (l/m r Q 2 R 4 ) 

advance ratio 
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<?S 

Vi 

% 

Ps 

Poo 

n 

( ) 4 
( )' 
(') 

( r 

NEL 

NOM 

rpm 


rotor solidity 

azimuthal location of the z-th blade 

z-th flap mode 

blade mass density (kg/m 3 ) 

air density (kg/m 3 ) 

rotor angular speed, (rad/sec) 

aerodynamic component of ( ) 

inertial component of ( ) 

time derivative of ( ) 

spatial derivative of ( ) 

approximately equal to 

number of aerodynamic elements 

number of modes each in flap bending, lag bending and torsion 
revolutions per minute 


1 Introduction 

In the early eighties, McNulty experimentally investigated the stability of isolated hin- 
geless rotors in forward flight under untrimmed conditions [1], The strengths of the ex- 
periment were structural simplicity and aerodynamically demanding conditions. In fact, the 
hub-flexure-blade assembly was intentionally designed to approximate closely root-restrained 
rigid flap-lag blades, and the database included dynamically stalled conditions with advance 
ratio as high as 0.55 and shaft angle as high as 20°. Many investigations right until to- 
day have been driven by the database in Ref. [I], and these strengths helped isolate the 
aerodynamic aspects of the stability problem in forward flight (Refs. 2—8), particularly the 
nonlinear dynamic stall aspects. However, the weaknesses of the experiment were low tip 
speed and the lack of a swash plate. For example, the collective pitch and shaft angle were 
the only control inputs, which were manually set prior to each run. This resulted in the 
unusual flight conditions of unrestricted tip path planes with very large negative coning an- 
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gles. To remove these weaknesses Maier et al. have been conducting experiments since the 
early ’90s on isolated hingeless rotors with torsionally soft blades [9]. Moreover, they also 
recognize the increasing need for a database on torsionally soft blades in forward flight [9]. 
Thus, stated explicitly, the primary objective of their tests was to generate trim and aeroe- 
lastic stability data for torsionally soft rotors operating at realistic rotational speeds. The 
rotor is operated under trimmed conditions of minimized 1/rev root flap moment (details to 
follow?). Thus, the database of Ref. [9] includes test results on lag-regressive mode damping, 
longitudinal and lateral cyclic pitch controls and steady root flap moment. 

With this background, we now spell out the objectives of this report: 

1. Develop an elastic flap-lag-torsion analysis for required trim conditions of minimized 
1/rev root flap moment. 

2. Predict trim results of longitudinal and lateral cyclic pitch controls, and steady root 
flap moment as well as stability results of modal frequencies and lag damping levels. 

3. Correlate these predictions with the database of Refs. [9] and [10] and thereby isolate 
the effects of quasisteady stall, dynamic stall and dynamic wake on trim and stability 
results. 

2 Experimental Rotor 

To facilitate discussion and presentation of results, we briefly describe the experi- 
mental rotor of Maier et al. [9]. It is a soft-inplane hingeless rotor with torsionally soft 
blades. The 7.5 -ft rotor has four blades of an NACA 0012 airfoil section. The blade chord 
is 3.4-in and has a rectangular plan form as shown in Fig. 1. Furthermore, the blades have 
two-degree precone, and zero-degree pretwist and droop. The rotor also has a provision to 
vary the blade precone. The blade mass-center, tensile, aerodynamic and elastic axes are 
nearly coincident with the control axis, which is at the quarter-chord point. Table 1 summa- 
rizes the rotor properties of the experimental model. The test-model blade including other 
details of the hub attachment are given in Figs. 1 and 2. Basically, the blade comprises four 
distinct regions. The first region is a hub section with very high stiffness values. The second 
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Figure 1: Hub-Flexure-Blade Assembly 

region is a root-flexure section, which accommodates the blade flap and lead-lag motions. 
The third region is a short transition section, which is relatively stiff and provides transition 
from the blade root-flexure to the airfoil section. The fourth section is the NACA 0012 
blade portion. Table 2 details the stiffness and mass distributions in each of these sections. 
The model is designed so that the test data correspond to stability of isolated rotors. The 
separation between the lag-regressive mode frequency and the lowest test-stand frequency is 
7.1 Hz, which is far above the lag regressive-mode frequency. 

The collective pitch angle (0 O ) and shaft tilt angle (or s ) are set prior to each run and 
are known parameters. The rotor is operated trimmed with lateral and longitudinal cj clic 
pitch controls. The cyclic pitch controls are exercised through a conventional swash plate 
mechanism that controls pitch on the blade root cuff at a location inboard of the flap-lag 
flexure motions (Fig. 1). According to Ref. [9], such an arrangement minimizes root-end- 
kinematic coupling effects and thus avoids any contamination of the aeroelastic stability data 
with kinematic nonlinearities. The rotor is operated at 1700 rpm, which in hover gives a 
Reynolds number of 1.2 x 10*^ and Mach number of 0.6 at the blade tip. 

In hover tests, the rotor is set at the desired speed and collective pitch angle. The 
shaft angle is set to zero degrees. Then, low-amplitude cyclic excitation is applied and 
the frequency of the excitation is adjusted until the lag response of the first blade reaches 
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Table 1: Details of the Experimental Rotor 


Number of blades 

4 

Airfoil section 

NACA 0012 

Hover blade-tip Mach number at 1700 rpm 

0.6 

Hover blade-tip Reynolds number at 1700 rpm 

1.2xl0 6 

Rotor radius, ft 

3.75 

Blade chord, in 

3.4 

Nonrotating fundamental flap frequency, Hz 

4.499 

Nonrotating fundamental lead-lag frequency, Hz 

14.405 

Nonrotating fundamental torsion frequency, Hz 

64.362 

Blade precone, deg 

0.0 and 2.0 

Blade pretwist, deg 

0.0 

Blade droop, deg 

0.0 

Blade sweep, deg 

0.0 
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Figure 2: Schematic of Hub-Flexure-Blade Assembly 

a maximum. To obtain a strong decaying transient signal, the magnitude of the cyclic 
excitation is adjusted till the maximum allowable lag response is reached; the maximum 
allowable is based on the blade structural limit. The excitation is stopped and the ensuing 
transient signal is recorded for 2 seconds. 

In forward-flight tests, the shaft angle is set first and the rotor speed is brought up to the 
desired value. Then, the collective pitch is set, wind tunnel air speed is increased slowly to 
the desired forward-flight speed and the cyclic pitch controls are adjusted to maintain low- 
oscillatory flapping loads. When the desired forward-flight speed is reached, the collective 
pitch is readjusted to get the desired value and the cyclic pitch controls are further adjusted 
to minimize the 1/rev flapping moment at the 12% radial station. Then, low-amplitude 
cyclic excitation is applied. The frequency and magnitude of the excitation are adjusted 
so that the maximum lead-lag response (but below the structural limit) is attained. The 
excitation is shut off and the ensuing transient is recorded for 2 seconds. The recorded 
signals are analyzed using the moving-block analysis technique to obtain modal damping 
and frequency. Reference [9] provides additional details on the model rotor design and test 
procedures. 
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3 Analysis 

3.1 Elastic Flap-Lag- Torsion Equations 

The flap-lag-torsion equations of motion are nonlinear partial differential equations, 
which are given in Ref. [6]. We use a Galerkin-type scheme to transform these partial 
differential equations into a set of ordinary differential equations in terms of generalized co- 
ordinates. Orthogonal, nonrotating normal modes are used, which are developed numerically 
with mass and stiffness distributions identical to those of the experimental rotor; see Table 
2. In-vacuo conditions are assumed; that is, without aerodynamics and the frequencies are 
obtained from an eigenanalysis. A Mykelstad-type approach is followed in generating the 
modes. As shown in Table 2, the stiffness distributions for the transition region (region 3) 
are not given because these distributions change continuously over this region and it was not 
possible to measure them [9]. In the present correlation work, we assume a linear variation 
for stiffness properties in that region while computing the mode shapes. Figures 3—5 show 
the nonrotating mode shapes for the first five flap bending, lag bending and torsion deflection 
modes. These modes are normalized with a tip deflection of one. The corresponding five 
nonrotating frequencies are tabulated in Table 3. Using these normal modes, the equations of 
motion are transformed into a set of modal equations in terms of the generalized coordinates. 
The Galerkin-type integrals associated with this transformation for spatial dependency are 
evaluated numerically and are given in Ref. [6]. 

3.2 Aerodynamics 

The airfoil aerodynamics is based on the quasisteady stall, dynamic stall and dynamic 
stall and wake theories, which are described in Ref. [6]. As in the experimental rotor, 
the analytical model has four blades of an NACA 0012 airfoil section. The airfoil-section 
quasisteady stall characteristics are determined following Ref. [11] at a Reynolds number 
of 1.8xl0 6 . For the lift model, the quasisteady lift coefficient C z , and extrapolated lift 
coefficient C z , are given by 

C Zl = a sin a cos a 0° < a < 180° U a ) 

C z ,=C Zi 0°<a<a° ss (lb) 
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Figure 3: Nonrotating Flap Bending Modes for the Experimental Rotor Blade 
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Figure 5: Nonrotating Torsion Modes for the Experimental Rotor Blade 


















C Za = a sin a 5S cosa ss a° s < a < 45° (lc) 

C Zi — a sin a ss cos a ss sin 2 q 45° < a < 135° (Id) 

C Zs = — a sin a s3 cos a ss 135° < a < 166° (ie) 

C, 3 =C Zl 166° < a < 180° (If) 


where a = 6.28 and = 14°. Figure 6 shows the quasisteady lift coefficient for 0° < a < 
360°. It also includes a comparison with the static test data from Ref. [11] and with the 
extrapolated linear-lift coefficient C Zl according to Eq. (la). 

For the drag model, the quasisteady drag coefficient Cd, and the constant drag coefficient 
C do are: 

C do = 0.01 ( 2a ) 

Cd 3 = 1 05 - (1.05 -Cd 0 ) cos 2a 0° < a < 360° (2b) 

The variations of these coefficients for 0° < a < 360° are shown in Fig. 7, which also 
includes test data from Ref. [11]. 

For the pitching moment model, the quasisteady moment coefficient C m , is given by 


C m , = C mo 0° < a < a° s (3a) 

C m , = C mo - 0.0582257tan~ l (a-a ss ) a° ss < a < 20° (3b) 

C ms = C mo - 0.55 sin(a - 20°) - 0.0842201 cos(a - 20°) 20° < a < 101.2941° (3c) 

C m s = C mo - 0.55641 cos[0.75(a - 101.2941°)] 101.2941° < a < 170° (3d) 

Cm, = 0.3461497[0.1(a- 170°)] 170° < a < 180° (3e) 


where C mo = 0.0. Figure 8 shows the variation of C m , from Eqs. (3a)— (3e) and the test 
data from Ref. [11]. 

The dynamic stall parameters used in the ONERA stall model are the same as those 
for an NACA 23012 airfoil section; these are given in the Appendix; for additional details, 
see Ref. [6]. To include wake effects, we use the finite-state three-dimensional wake model, 
which is based on six harmonics and three radial shape functions for each harmonic; this 
leads to 39 wake states. Although a detailed convergence study with respect to the number 
of harmonics and radial shape functions has not been carried out in forward flight, the wake 
model should provide results close to the converged values for Cj/^s < 0.05. 
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Figure 6: Variation of Airfoil Lift Coefficient for an NACA 0012 Airfoil Section 
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Figure 8: 




3.3 Trim Analysis 

As detailed earlier, the collective pitch and shaft tilt angles are known control inputs. 
Therefore, trim analysis per se refers to finding the lateral and longitudinal cyclic control 
inputs for minimized 1/rev flap moment at the 12% radial station and to finding the cor- 
responding initial conditions for the periodic responses. The trim analysis is based on the 
shooting method with damped Newton Iteration, which generates the FTM as a byproduct. 
The modal damping levels and frequencies are obtained from the eigenvalues and eigenvec- 
tors of the FTM [12]. This analysis still requires two additional trim equations that satisfy 
the required trim conditions of minimized 1/rev flap moment at 0.12/?. These trim equations 
are given in the next section. 


3.3.1 Flap Moment Equations 


Consider a generic material point on the blade at a radial location x from the hinge 
offset. As shown in Fig. 9, the generic point is subjected to inertial, aerodynamic and 
centrifugal forces. It is expedient to express the components of each of these forces in the 
undeformed blade coordinate system xyz. Let £ u and £ w represent the total horizontal 
and vertical forces parallel and perpendicular to the undeformed x-coordinate, respectively. 
Similarly, let M v represent the total moment acting parallel to the undeformed y-coordinate 
and let ej represent the 0.12 R radial location about which the flap moment is computed. 
Following Ref. [13], we predict the flap moment by the force-integration approach, in which 
the sectional forces and moments are integrated over the blade span to obtain the total 
moment. Therefore, the total flap moment is given by 


M(e f ,t)= f {£ u [w(x) - w(ef)} - £ w (x - e f ) + M u } dx (4) 

Je, 

Now, we express the forces and moment in Eq. (4) as a sum of inertial and aerodynamic 
components. The inertial component includes contribution from the centrifugal forces as 


well. Thus we have 


£ u = Cl + £t A 


= £l + C A 


w ’ 


M v = M[ + 


(5) 
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where the superscripts I and A indicate that the contributions are of inertial and aerody- 
namic origin. Substituting Eq. (5) in Eq. (4), we get 

M{e f ,t) = f [C[ [u7(x) - w{e f )] - C l w {x -e f ) + M[) dx 

Je f 

Je f 

— yW ( Cf ) inertial H" ) aerodynamic ( 6 ) 

The expressions for sectional inertial forces, and C l w , can be derived by applying 
Newton’s second law: 


C 1 = C[ i + C [ j + Cl k = JJ Ps a dr) d£ (7) 

Substituting for the blade acceleration relative to an inertial frame and performing the cross- 
sectional integrals, we get [6] 

C ! u = —m J— 2vQ — (z + e h )n 2 + /? pc u>fl 2 j (8) 

Cl = —m To + 2/? pc flv + /3 pc £l 2 (x + ) j (9) 

In deriving the above equations, we neglect the axial deformation. 

Similarly, the sectional inertial moments can be derived from 

M 1 = Ml i + Ml j + Ml k = - Jj p s s x adrj d£ (10) 


where s represents the moment arm (see Ref. [6]). Sustituting for s and a and performing 
the cross-sectional integrals, we get 

Ml = -m [-li' + flV] {(k 2 m3 - fc^Jsinflcosf} 

— m + 2fl0 / + + fl 2 /? pc ] jk m2 sin 2 0 + cos 2 # j (11) 

Substituting Eqs. (8), (9) and (11) for the inertial part in Eq. (6) and after nondimension- 
alization, we get 

M !^ = AT(e,,^) 1 «„ w =- /‘{2m >,[»(*)- 

p^biRR 4 7 Je f 
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( 12 ) 


+ m(x + e h )[w(x) - w(e f )\ 

— 0 pc m w [u>(:r) - w(ej)\ + m w(x - e/) 4- 2/3 pc m v(x — e/) 

+ / 3 pc m(x — e/)(x + e&) + mk^v' s\n 9 cos 9 — mfc^ t/sin 9 cos 9 
+ mk^w sin 2 9 — mk^w 1 sin 2 9 — 26' mk^ sin 2 9 

- mk^flpcs'm 2 9} dx 


Now, substitution of modal expansions for w and v in Eq. (12) gives 


6a 


NOM NOM 


ai = -{£ £ [2 FMlyW-PpcFMiWW 

1 »=1 j= 1 


NOM 

+ £ [FM 3 W t + FMfW l +2l3p C FM^V l 

4=1 

NOM 

+ £ [F Mf sin 9 cos 9Vi — FMf sin 9 cos 9Vi 

»= i 
NOM 

+ Y [p sin 2 9W{ — F Mj sin 2 0VE,J 

4 = 1 

+ fi pc FM 8 - (29 + p pc ) sin 2 9F M 9 } 


(13) 


The aerodynamic contribution to flap moment is obtained from sectional aerodynamic 
forces and moments, which are calculated at the mid point of each aerodynamic element. 
Hence, we use a force-summation scheme instead of integrating the forces and moments along 
the length of the blade. Thus, in terms of sectional aerodynamic forces and moments, the 
flap moment due to aerodynamics can be obtained as 


J^d (ey , 0 ) aerodynamic 

Poc bWR* 


»4d ( , 0 ) aerodynamic 


iVC/Ij 

Y [^u { w ( x j) ~ ™( e /)} - £i x i + yV M n 


(14) 


In the above equation, x f n and Ax n represent the radial distance of the mid point from the 
0.12/? radial station and the length of the n-th aerodynamic element, respectively. Expressing 
Eq. (14) in terms of sectional aerodynamic forces and moments, we get 


C* = — (v'L v + w'L w ) 


(15) 
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( 16 ) 



M 


A 

V 


v'M t 


(17) 


where L v , L w and are the sectional aerodynamic forces and moments described in Ref. 
[6]. Substitution of Eqs. (15)— (17) into Eq. (14) gives 

NEL 

Ad (c / 1 0) aerodynamic )3 [ y 10 Lyj) [tu(xy ) 10^6 j )] 

n=l 


- L w x j + v'M<f,] n Ax n 


(18) 


Substituting modal expansions, we rewrite the above equation as 

NEL NOM NOM 

d„™k = EfE E l-FM!°V i W 1 L„„-FMVW,W i L, 

n = 1 1=1 j=l 

NOM 

+ ^ FM? 2 V t M K -L Wn x fn } 

1=1 


(19) 


Similarly, with the substitution of Eqs. (13) and (19) in Eq. (6), the total flap moment at 
the radial location e/ is given by 
c NOM NOM 

M(e,) = -{£ E [2FM! i VlW i -/) r 'FMlW i W,] 

7 i=i j= i 

NOM 


+ £ [FM 3 W> + FMfWi + 2^ c FMtV l \ 

i=i 

NOM 

+ J2 [ FM f sin 9 cos 8V t - FM? sin 0 cos 0V{] 


1 = 1 
/VOAf 


+ £ [FM, 7 sin 2 9W, - FM? sin 2 8W t 

t=i 

+ fi pc FM % - (28 + fi pc ) sin 2 8FM 9 } 

NEL NOM NOM 

+ E(E Z[-FM$V i W i L..-FM$W i W i L, 

n=l i=l j=l 
jVOM 

+ E FMPViMt, - ( 

1 = 1 

The Galer kin-type integrals accounting for spatial contributions are detailed in Ref. [6]. 
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3.3.2 Trim Equations 


Equation (20) gives the total flap moment acting at the radial station e/; that is, at the 
12% radial station. To obtain the 1/rev components of flap moment, we expand Eq. (20) in 
terms of Fourier series as 

00 

M(e f ,ip) = Mo + [Mnccos(nip) + M ns sin(ni/>)] 

n=l 

(21) 

For n — 1, we get 


M (e /? ip) = M 0 {ef) + Mi c {e f ) cos ip + M u {ef) sin ip 

(22) 

where 


1 [2* 

Mic{ej) = - M(e fl ip) cos ip dip 

7 T Jo 

(23) 

1 /2 7T 

M n(e f ) = - / M(e f ,ip) smipdip 

7T J 0 

(24) 

Therefore, the required trim conditions are 

l [2ir 

Mic(ef) = ~ 1 M(e/, ip) cos ip dip = 0 

7T JO 

(25) 

Mi s (e/) = - / M{e f ,ip) sin ip dip = 0 

(26) 


7 r Jo 


Thus, Eqs. (25) and (26) represent the trim equations, which are solved together with the 
response periodicity conditions. 


3.3.3 Steady Flap Moment Based on the Curvature Method 

Concomitant to the force-integration method, the steady flapping moment at the radial 
location e f is also estimated by the curvature method. That is, 


-7-7/ , x _ 

,M( e/ ,0) n2R 2 


d 2 


w 


dx 1 


In terms of mode shapes and generalized coordinates, Eq. (27) is expressed as 

NOM 

M(ef,4>)= E 

l-l 

Therefore, the steady flap moment is given by 

M 0 (ej) = 77 - / M(e f ,ip) dip 

Ztt Jo 


(27) 


(28) 


(29) 
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4 Correlations 


We begin the correlation with the nonrotating and rotating frequencies of the blade 
modes in flap bending, lag bending and torsion. We then correlate with the lag regressive- 
mode damping, control settings of lateral and longitudinal cyclic pitch angles and finally root 
flapping moment at the 12% radial station. The baseline parameters for these predictions 
are identified in Table 4. 

Table 3 shows a comparison between the data and predicted nonrotating frequencies for 
flap bending, lag bending and torsion modes; the data are available for three fundamental 
flap-lag-torsion modes as well as for the second and third modes in flap bending. The 
correlation is very good in that the maximum error in the prediction with respect to the 
data hardly exceeds 5.5%. Figure 10 shows the variation of the rotating frequencies with 
rotor rotational speed f l for the first two flap-lag-torsion modes. For example, in forward 
flight 0 = 1700 rpm and the corresponding predictions of the dimensionless (dimensional) 
rotating frequencies for the fundamental flap-lag-torsion modes are 1.13 (32.02 Hz), 0.7 
(19.83 Hz) and 2.86 (81.03 Hz), respectively. The rotating-frequency data are available 
over a range of fl values for the first lag mode only, and here as well, the correlation is 
very good. For completeness, Fig. 10 also includes the nonrotating-frequency data on the 
fundamental flap-lag-torsion and second-flap bending modes. 

Compared to the frequency data in Table 3 and Fig. 10, the data on damping level, 
cyclic pitch controls and root flapping moment are far more extensive [9, 10]. Accordingly, 
we present the correlation with these data in three phases: lag regressive-mode damping in 
the first phase (Figs. 11-14), lateral and longitudinal cyclic pitch controls in the second 
phase (Figs. 15-22), and finally root flapping moment in the third phase (Figs. 23-30). 
Following Maier [10], we select four typical test cases with a two-degree precone, for which 
the updated database is available: three-degree collective with three cases of shaft angle 
a s = 0°, 3° and 6°, and 5.9-degree collective with shaft angle a s = 6°. The analytical 
predictions are from the quasisteady stall, dynamic stall, and dynamic stall and wake aero- 
dynamic theories. As seen from the data in Fig. 11 (0o = 3° and a s = 0 ), the damping 
increases slowly up to an advance ratio of 0.2, and thereafter it more or less decreases with 
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Table 4: Baseline Parameters for the Experimental Rotor 


Lift curve slope, a (1 /rad) 6.283 

Lock number, 7 (based on a = 6.283) 7.5896 

Profile drag coefficient, Cd 0 0.01 

Camber of NACA 0012 airfoil section, a c 0.0° 

Lift coefficient at a = 0, Ci 0 0.0 

Pitching moment coefficient at a = 0, C ma 0.0 

Quasisteady stall angle of NACA 0012 airfoil section, a ss ±14° 
Rotor solidity, a s 0.0962 

Hinge offset, th (in) 4.7 

Blade root cutout, r c (in) 13.76 

Blade precone, /3 pc (deg) 2 

Number of aerodynamic elements, N EL 5 

Number of modes, NOM 5 


increasing fi. Thus, the data exhibit a convex trend of damping variation with fi. Overall, 
this trend of the data is well predicted by all three theories. Furthermore, the thrust level 
Ct/vs is low throughout the data range (0 < fi < 0.31); in fact, it is less than 0.03, and for 
fi > 0.2, it slowly decreases with increasing advance ratio. Given the low-thrust condition, 
this is expected. The comparative aspects of the correlation by three theories are instructive 
as well. For example, despite some overprediction for fi < 0.05, the dynamic stall theory 
provides good correlation throughout. Similarly, the quasisteady stall theory closely follows 
and underpredicts the dynamic stall theory. This underprediction, though marginal, is con- 
sistently observed throughout the data range. Compared to these two theories, the dynamic 
stall and wake theory provides overall better correlation in that it picks up the finer varia- 
tions of the experimental data. For example, for 0 < fi < 0.05, it predicts a nearly constant 
damping level and agrees with the data better, and for ft > 0.15, it gradually follows the 
dynamic stall theory. The minor differences between the quasisteady stall and dynamic stall 
theories are due to unsteady lift effects; differences between dynamic stall and dynamic stall 
and wake theories are due to dynamic wake effects. 

In Fig. 12, #0 = 3° and a s = 3°. Here as well, the data show a convex trend of 
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the damping level. The thrust level, which is about 0.01 in hover, increases to 0.022 for 
/r « 0.1 and thereafter decreases to 0.013 at fi = 0.31. Since the thrust level is low, all three 
aerodynamic theories provide good correlation. In particular, the dynamic stall and wake 
theory reduces the overpredictions of the quasisteady stall and dynamic stall theories for 
H < 0.1 and gradually follows the other two theories for fi > 0.2. 

Now, we come to Fig. 13 for 6 0 — 3° and cr s = 6°. As seen from the data, the damping 
level very slowly increases with increasing advance ratio for fi < 0.1 and thereafter, it 
sharply decreases. The thrust level is about 0.01 in hover and increases to a maximum value 
of 0.017 at /z = 0.15 and thereafter decreases with increasing fi ; in fact, it is nearly zero at 
H = 0.31. The predictions from the quasisteady stall and dynamic stall theories show slowly 
increasing damping up to /z « 0.1 and thereafter decreasing damping. However, the rate of 
decrease in damping is small. Moreover, these two theories somewhat overpredict damping 
for 0 < fi < 0.05. By comparison, the dynamic stall and wake theory improves the correlation 
for n < 0.15, and for /z > 0.15, it slowly merges with the other two theories. This merging 
is expected since the thrust level, which is low to begin with, decreases with increasing fi for 
H > 0.15 and more or less reduces to zero at /z = 0.31. In summary, the correlation from the 
dynamic stall and wake theory is adequate for fi < 0.2 and needs quantitative corrections 
for fi > 0.2, which merits further investigation. 

Figure 14 is for 6 0 = 5.9° and a s = 6°, and the data are available from a near-hovering 
condition with fi = 0.04 to a fairly high-speed condition with fi = 0.36. Despite some data 
scattering around ft = 0.05 and 0.35, the data show that the lag damping increases for 
0.04 < fi < 0.15 and decreases for fi > 0.15. We also mention that the thrust level increases 
from 0.034 in hover to 0.045 at /z « 0.1 and decreases to 0.02 at fi = 0.36. Overall, all three 
theories predict the trend of the data. However, the dynamic stall theory and the quasisteady 
stall theory (except at fi fa 0.15) overpredict damping consistently. The differences between 
the quasisteady stall and dynamic stall theories are mainly due to unsteady lift effects, which 
increase with increasing collective pitch; see Fig. 13. The dynamic stall and wake theory 
significantly reduces the overpredictions of the quasisteady stall and dynamic stall theories 
and thereby improves the correlation. Overall, the dynamic stall and wake theory provides 
good correlation. 
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In Figs. 15-17, we present the correlation for lateral cyclic-pitch angle 9 C for o? 3 = 0°, 3° 
and a s = 6°, respectively. In all three figures, the data show known trends: 9 C is zero at 
/i = 0.0, increases suddenly around fi . = 0.05 and thereafter slowly decreases with increasing 
advance ratio. The rate of decrease is higher in Fig. 17 than in Figs. 15 and 16. The 
predictions from the quasisteady stall and dynamic stall theories are nearly the same for 
all three values of a s . And these predictions show that the lateral cyclic pitch required to 
minimize 1/rev flapping moment is negative and that its magnitude monotonically increases 
with increasing advance ratio. Although these trends of the predictions are at best consistent 
with those of the data for fi > 0.05, overall they are quantitatively inaccurate. But inclusion 
of wake effects dramatically improves the correlation, both quantitatively and qualitatively. 
In particular, the dynamic stall and wake theory predicts the sudden increase around [i = 
0.05; that is, in the transition regime when the flow over the rotor disk' is associated with a 
large amount of shed and trailing vorticities. With increasing fi, the dynamic stall and wake 
theory predicts the decreasing trend of the data. However, this rate of decrease is relatively 
higher and, hence, the correlation needs further quantitative correction for // > 0.2. 

In Fig. 18, 6 0 = 5.9° and a s =6°. Here as well, the data show a variation for 0 C similar to 
that in Figs. 15—17; that is, suddenly increasing at = 0.05 and thereafter decreasing with 
increasing advance ratio. As seen from Fig. 18, both the dynamic stall and quasisteady stall 
theories fail to predict this variation and are not acceptable. By comparison, the dynamic 
stall and wake theory not only predicts this variation but also provides good correlation 
throughout. In summary, Figs. 15-18 show that it is important to include wake effects in 
the prediction of lateral cyclic pitch angle 0 C and that the correlation in Figs. 16 and 17 
merits further improvement at high speed conditions, say, ji > 0.25. 

Now, we show the correlation for longitudinal cyclic pitch angle 9 S for four cases, three- 
degree collective with a 3 = 0°, 3° and 6° in Figs. 19-21, respectively, and 5.9°-collective 
with cx s — 6° in Fig. 22. The data are available from hover to an advance ratio of 0.31 
and show that negative 9 S is required for the present trim condition of minimizing 1/rev 
root flapping moment. As for the details, we begin with Fig. 19; the data show that 9 S 
is nearly zero up to n « 0.05 and that thereafter it increases (negative) with increasing p. 
The predictions from the quasisteady stall and dynamic stall theories show that 9 S essentially 
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increases (negative) with increasing fi, and the predictions from these two theories are nearly 
identical. Compared to the data, however, this rate of increase is milder. Moreover, the 
theories fail to capture the finer details in the variations of 9 S with fi at low advance ratios, 
say for /t < 0.1. By comparison, the dynamic stall and wake theory predicts that 0 3 is nearly 
constant for fi < 0.075. Thereafter, however, it basically follows the other two theories and, 
all in all, requires a bit of quantitative improvement for fi > 0.1. Similar trends are observed 
for o s = 3° in Fig. 20. In summary, as seen from Figs. 19 and 20, the quasisteady stall 
and dynamic stall theories provide satisfactory correlation and the dynamic stall and wake 
theory follows the other two theories with noteworthy improvements for ft < 0.05 and some 
overpredictions for 0.05 < fi < 0.1. 

In Fig. 21 the data show a similar variation for 9 S as in Figs. 19 and 20. Here as well, the 
predictions from the quasisteady stall and dynamic stall theories agree, and they correlate 
well with the data. The dynamic stall and wake theory brings the theory closer to the data 
at ft = 0.05 but it also takes the theory somewhat away from the data at fi = 0.1. For 
H > 0.1 it virtually merges with the other two theories. Overall, all three theories provide 
good correlation. 

The data in Fig. 22 show that 9 S (negative) increases with increasing advance ratio. 
The predictions from the quasisteady stall and dynamic stall theories are nearly identical 
and provide adequate correlation. The dynamic stall and wake theory slightly improves 
the correlation for ft > 0.25, but compared to the other two stall theories, it also slightly 
overpredicts for fi < 0.25. Nevertheless, it also provides adequate correlation overall. 

The next four figures show the correlation for the root flap moment from the curvature 
method. The three-degree collective case with a s = 0°, 3° and 6° is covered in Figs. 23-25, 
respectively, and the 5.9-degree collective case with a s = 6° is covered in Fig. 26. The 
data show that the flap moment increases with increasing fi up to a certain value of fi, say, 
around 0.2 in Figs. 23 and 24 and around 0.15 in Figs. 25 and 26, and that it decreases 
thereafter; that is, a convex-type variation. Surprisingly, the predictions from all three 
theories are nearly identical. Although these three theories predict the trend of the data, 
they underpredict (negative) flapping moment throughout the data range except in Fig. 25 
(i 9 0 = 3° and a s = 6°) for fi > 0.25 and in Fig. 26 at p ~ 0.36. 
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Figures 27 — 30 are based on the force integration method. Here as well, the predictions 
from the three aerodynamic theories are nearly identical and they consistently underpredict 
the flapping moment for all values of fi. Nevertheless, they predict the trend of the data for 
all four cases: a s = 0°, 3° and 6° with 8 0 = 3° and a a = 6° with 8 0 = 5.9°. Overall, the 
required quantitative improvements in Figs. 23—30 merit further investigation. 


5 Toward Improving the Correlations 

We are pursuing a two-pronged approach toward improving the preceding correlation. 
The first approach is based on improving the dynamic stall characteristics and the results are 
presented in the next section. The second approach is based on adapting a transient-response 
analysis [14] somewhat similar to the moving block analysis. It is still in a developmental 
stage and the research is continuing. 

5.1 NACA0012 Dynamic Stall Characteristics 

Thus far, the correlations in Figs. 11-30 are based on the dynamic stall parameters of 
an NACA 23012 airfoil section used in Ref. [3]. It needs to be emphasized that the database 

on dynamic stall parameters is extremely limited and that the state of the art does not 

permit an accurate description of these parameters, no matter which airfoil. (Throughout, 
the quasisteady stall parameters of an NACA 0012 airfoil section are used.) Given this 
background, it is instructive to study the sensitivity of the preceding correlations to the 
dynamic stall parameters. That study is presented in Figs. 31—50, in which the correlations 
are based on the dynamic stall parameters of an NACA 0012 airfoil section used in Ref. 
[8]. These correlations are also presented in three phases: Figs. 31—34 on lag-regressive 
mode damping levels, Figs. 35-38 on lateral cyclic pitch control input, Figs. 39-42 on 
longitudinal cyclic pitch control input, and Figs. 43—50 on root flap moment. 

Figures 31-33 show the sensitivity of lag-regressive mode damping to dynamic stall 

parameters for d 0 = 3° with a., = 0°, 3° and 6°, respectively, and Fig. 34 for B Q = 5.9° 

and <* s = 6°. Overall, Figs. 31-34 show that the sensitivity of the damping predictions 
to changes in the dynamic stall parameters is negligible. The remaining sets of correlations 
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further demonstrate this negligible sensitivity, Figs. 35—42 with respect to the lateral cyclic 
pitch control and longitudinal cyclic pitch control angles and Figs. 43—50 with respect to 
the root flap moment. 

6 Concluding Remarks 

The preceding correlation covers a comprehensive database (0.0 < fj, < 0.36, 0° < 9 S < 
5.9°, 0° < a s < 6°) on lag damping, lateral and longitudinal cyclic pitch angles and root 
flap moment. Overall, the dynamic stall and wake theory provides fairly good correlation. 
A major finding is that inclusion of dynamic wake dramatically improves the correlation for 
the lateral cyclic pitch angle 0 C . We are continuing this correlation effort to cover additional 
data received recently from the AFDD and to adapt a transient-response-analysis method 
of Ref. [14], which provides a means of predicting damping levels and frequencies without 
perturbations as required in Floquet analysis. 
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Figure 10: Fan Diagram for Rotating Frequencies of the Experimental Rotor 
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Figure 1 1: Effects of Aerodynamic Modeling on Lag Damping Correlation 

for 0 Q = 3° and a s = 0° 
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Figure 12: Effects of Aerodynamic Modeling on Lag Damping Correlation 

for 0 O = 3° and a s = 3° 
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Figure 13: Effects of Aerodynamic Modeling on Lag Damping Correlation 

for 0„ = 3° and a = 6° 
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Figure 14: Effects of Aerodynamic Modeling on Lag Damping Correlation 

for 0„ - 5.9° and a = 6° 
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Advance Ratio (}x) 

Figure 15: Effects of Aerodynamic Modeling on Lateral-Cyclic Pitch Correlation 

for 0 O = 3° and a s = 0° 
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Figure 16: Effects of Aerodynamic Modeling on Lateral-Cyclic Pitch Correlation 

for 0 () = 3° and a s = 3° 
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Figure 17: Effects of Aerodynamic Modeling on Lateral-Cyclic Pitch Correlation 

for 0 () = 3° and a s - 6° 
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Advance Ratio (|i) 

Figure 18: Effects of Aerodynamic Modeling on Lateral-Cyclic Pitch Correlation 

for 0 O = 5.9° and a s = 6° 
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Figure 19: Effects of Aerodynamic Modeling on Longitudinal-Cyclic Pitch 
Correlation for 0„ = 3° and a = 0° 
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Figure 20: Effects of Aerodynamic Modeling on Longitudinal-Cyclic Pitch 
Correlation for 0 O = 3° and a s = 3° 
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Figure 21: Effects of Aerodynamic Modeling on Longitudinal-Cyclic Pitch 
Correlation for 0 O = 3° and a s = 6° 
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Advance Ratio (|i) 

Figure 22: Effects of Aerodynamic Modeling on Longitudinal-Cyclic Pitch 
Correlation for 0 O = 5.9° and a s = 6° 
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Figure 23: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Curvature Method for 9 0 = 3° and a s = 0° 
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Figure 24: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Curvature Method for 0 O = 3° and a s = 3° 
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Figure 25: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Curvature Method for 0 Q = 3° and a s = 6° 
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Figure 26: Effects of Aerodynamic Modeling on Flapping Moment Correlaion 
from the Curvature Method for 0„ = 5.9° and a = 6° 
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Figure 27: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Force-Integration Method for 0 O = 3° and a s = 0° 
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Figure 28: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Force-Integration Method for 0 O = 3° and a s = 3° 
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Figure 29: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Force-Integration Method for 0 () = 3° and a s = 6° 
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Figure 30: Effects of Aerodynamic Modeling on Flapping Moment Correlation 
from the Force-Integration Method for 0 O = 5.9° and a s = 6° 
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Figure 3 1 : Effects of Dynamic Stall Parameters on Lag Damping from 
the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 0° 
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Figure 32: Effects of Dynamic Stall Parameters on Lag Damping from 
the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 3° 
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Figure 33: Effects of Dynamic Stall Parameters on Lag Damping 
from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 6 ( 
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Figure 34: Effects of Dynamic Stall Parameters on Lag Damping from the 
Dynamic Stall and Wake Theory for 0 () = 5.9° and a = 6° 
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Figure 35: Effects of Dynamic Stall Parameters on Lateral-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 0° 
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Figure 36: Effects of Dynamic Stall Parameters on Lateral-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 0 O = 3° and oc. = 3° 
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Figure 37: Effects of Dynamic Stall Parameters on Lateral-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 9 0 = 3° and a s = 6° 
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Figure 38: Effects of Dynamic Stall Parameters on Lateral-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 0 O = 5.9° and a s = 6° 
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Figure 39: Effects of Dynamic Stall Parameters on Longitudinal-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 0° 
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Figure 40: Effects of Dynamic Stall Parameters on Longitudinal-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 0 () = 3° and a s = 3° 
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Figure 41: Effects of Dynamic Stall Parameters on Longitudinal-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theroy for 0 O = 3° and a s = 6° 
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Figure 42: Effects of Dynamic Stall Parameters on Longitudinal-Cyclic Pitch Angle 
from the Dynamic Stall and Wake Theory for 0 O = 5.9° and a s = 6° 
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Figure 43: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Curvature Method from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 0° 
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Figure 44: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Curvature Method from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 3 C 
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Figure 45: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Curvature Method from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 6° 
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Figure 46: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Curvature Method from the Dynamic Stall and Wake Theory for 0 O = 5.9° and a s = t 
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Figure 47: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Force-Integration Method from the Dynamic Stall and Wake Theory for 0 O = 3° and oc s = O' 
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Figure 48: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Force-Integration Method from the Dynamic Stall and Wake Theory for 0 Q = 3° and a s 
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Figure 49: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Force-Integration Method from the Dynamic Stall and Wake Theory for 0 O = 3° and a s = 6° 
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Figure 50: Effects of Dynamic Stall Parameters on Flapping Moment Based on the 
Force-Integration Method from the Dynamic Stall and Wake Theory for 0 O 5.9° and a s = 6 1 
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Appendix 

Dynamic Stall Parameters 
NACA 23012 Airfoil Section(Ref. 6) 


A = 0.2 

(30a) 

_ dC *< a 
da 2 

(30b) 

w = 0.10 + 0.023(|q°| — 13°)u(|a°| — 13°) if |a°| < 21.7° 


= 0.3 if |a°| > 21.7° 

(30c) 

, 0.105 

a — 

w 

(30d) 

e = 2 — 5.1 tan -1 {l.21(|a 0 | — 13°)} x u(|a°| — 13°) 


where u is the unit step function 

(30e) 

a d = 0.32 

(31a) 

r d = 0.2 + 0.1(AC z ) 2 

(31b) 

E d = -0.015 (AC,) 2 

(31c) 

a m = 0.25 + 0.1 {AC z f 

(32a) 

r m = 0.2 + 0.2 (AC,) 2 

(32b) 

E m = 0.01 (ACJ 2 

(32c) 

NACA 0012 Airfoil Section (Ref. 8) 


il 

0 
>— 1 

01 

(33a) 

c dc, a 

da 2 

(33b) 
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w 2 (l+d 2 ) = [ 0 . 2 + 0.1 (AC',) 2 ] 2 

(33c) 

2dw = 0.25 + 0.1 (AC,) 2 

(33d) 

e = —0.6 (AC,) 2 

(33e) 

a d = 0.32 

(34a) 

r d = 0.2 + 0.1 (AC,) 2 

(34b) 

E d = -0.015 (AC,) 2 

(34c) 

o m = 0.25 + 0.1 (AC’,) 2 

(35a) 

r m = 0.2 + 0.1 (AC,) 2 

(35b) 

E m - -0.6 (A C z f 

(35c) 
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